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We calculate analytically the proportionality constant in 
the pressure law of a membrane between parallel walls from 
the strong-coupling limit of variational perturbation theory 
up to third order. Extrapolating the zeroth to third ap- 
proximations to infinity yields the pressure constant a = 
0.0797149. This result lies well within the error bounds 
of the most accurate available Monte-Carlo result a MC = 
0.0798 ± 0.0003. 



I. MEMBRANE BETWEEN WALLS 

The violent thermal out-of-plane fluctuations of a 
membrane between parallel walls generate a pressure p 
following the law 
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(1.1) 



whose form was first derived by Helfrich Q. Here, n 
denotes the elasticity constant of the membrane, and d 
the distance between the walls. The exact value of the 
prefactor a is unknown, but estimates have been derived 
from crude theoretical approximations by Helfrich JjJ and 
by Janke and Kleinert H which yielded 



c4 h w 0.0242, 



,th 
*JK 



0.0625. 



(1.2) 



More precise values were found from Monte-Carlo simu- 
lations by Janke and Kleinert Q and by Gompper and 
Kroll p] which gave 



MC 



0.079 ±0.002, a$£ « 0.0798 ±0.0003. (1.3) 



In a previous work Q, a systematic method was devel- 
oped for calculating a with any desired high accuracy. 
Basis for this method is the strong-coupling version of 
variational perturbation theory 0|. In that theory, the 
free energy of the membrane is expanded into a sum of 
connected loop diagrams, which is eventually taken to 
infinite coupling strength to account for the hard walls. 
As a first approximation, an infinite set of diagrams was 
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calculated, others were estimated by invoking a mathe- 
matical analogy with a similar one-dimensional system of 
a quantum mechanical particle between walls. The result 
of this procedure was a pressure constant 
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0.0771063.. 



(1.4) 



very close to (|l -3[) . 

It is the purpose of this paper to go beyond this esti- 
mate by calculating all diagrams up to four loops exactly. 
In this way, we improve the analytic approximation (1.4) 
and obtain a value 



a th w 0.0797149, 



(1.5) 



which is in excellent agreement with the precise MC value 
a™£ in Eq. (Oh. 



II. THE SMOOTH POTENTIAL MODEL OF THE 
MEMBRANE BETWEEN WALLS 

To set up the theory, we let the membrane lie in the 
x-plane and fluctuate in the z-direction with vertical dis- 
placements y(x). The walls at z = ±d/2 restrict the 
displacements to the interval ip 6 (—d/2,d/2). Near 
zero temperature, the thermal fluctuations are small, 
¥>(x) ~ 0. The curvature energy E of the membrane 
has the harmonic approximation H 



E=^k I dx 2 [<9VW] 2 



(2.1) 



The thermodynamic partition function Z of the mem- 
brane is given by the sum over all Boltzmann factors of 
field configurations y(x) 
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This simple harmonic functional integral poses the prob- 
lem of dealing with a finite range of fluctuations. This 
problem is solved by the strong-coupling theory of Ref . Q| 
as follows. 

If the area of the membrane is denoted by A, the par- 
tition function (2.6) determines the free energy per area 
as 
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By differentiating / with respect to the distance d of the 
walls, we obtain the pressure p = —df/dd. 



From Refs. the harmonic free energy per unit area 

f m 2 is known as 
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(2.9) 



A. Smooth Potential Adapting Walls 

We introduce some smooth potential restricting the 
fluctuations y(x) to the interval (—d/2, d/2), for instance 

d 2 TT 
V(<p(x)) = m 4 — tan 2 -ip(x) 
d 
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= mV(x) + ^V fat (¥>(x)), 



(2.4) 



where we have split the potential into harmonic and in- 
teracting part 



^„t(^(x)) - m 4 [e 4 </(x) +e 6 V(x) + 
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with e 4 = 2/3, e 6 = 17/45, e 8 = 62/315, 
are left with the functional integral 



(2.5) 
Thus we 



XV (x) exp ( — 



d 2 x- 



{[dV(x)] 2 +mV(x) + 

JVtat^x))}), (2.6) 

where we have set k = fc^T = 1. After truncating the 
Taylor expansion around the origin, the periodicity of the 
trig ono metric function is lost and the integrals over (f(x) 
in ( |2.2| ) can be taken from — oo to +oo. The interact- 
ing part is treated perturbatively. Then, the harmonic 
part of V(ip(x)) leads to an exactly integrable partition 
function Z rn i . The mass parameter m is arbitrary at 
the moment, but will eventually taken to zero, in which 
case the potential V(ip(x)) describes two hard walls at 
ip = ±d/2. 

We shall now calculate a perturbation expansion for Z 
up to four loops. This will serve as a basis for the limit 
Tfi — > 0, which will require the strong-coupling theory of 
Rcf. §. 

B. Perturbation Expansion for Free Energy 

The perturbation expansion proceeds from the har- 
monic part of Eq. ( 2.20| ): 



z m 2 = <j> xy(x) « 



-A m2 [<p]/2 _ -Af 77 



(2.7) 



with 



d 2 i{[9VW] 2 +mV(x)}. (2 



The harmonic correlation functions associated with (2.7) 
are 

(Oi(^(xi))0 2 (<p(x 2 )) ■••)m 2 = 
' / 'l?< /3 (x)0 1 V(x 1 ))0 2 (^(x 2 ))-.. e ^M/2 (2 . 10 ) 



where the functions Oi((p(x.j)) may be arbitrary polyno- 
mials of ip (x 3 ). The basic harmonic correlation function 



G m 2(xi,x 2 ) = (<ys(xi)v?(x 2 )) r 



(2.11) 



deter mines, by Wick's rule, all correlation functions 
( |2.10|) as sums of products of ( 2.11 ): 



^ G m 2(xp(i),xp( 2 )) • • • G m 2(x P( 

), (2-12) 

pairs 

where the sum runs over all pair contractions, and P de- 
notes the associated index permutations. The harmonic 
correlation function (2.11) is in momentum space 



G m2 (k) 
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k 2 + im 2 



(2.13) 



thus being proportional to the difference of two ordi- 
nary correlation functions (p 2 — /i 2 ) -1 with an imaginary 
square mass /x 2 = ±im 2 . From their known x-space form 
we have immediately 

G m 2 (xi , x 2 ) = G m 2 (xi - x 2 ) 

K (Vim\xi - x 2 |) - K (y/^im\xi — x 2 |) , 

(2.14) 
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where Kq{z) is a modified Bessel function || Sec. 8.432]. 
At zero distance, the ordinary harmonic correlations 
are logarithmically divergent, but the difference is finite 
yielding G m 2(0) = l/8m 2 . 

We now expand the partition function (|2.6|) in powers 
of gVi n t((j)(x)) , where g = ir 2 /d 2 , and use (2.1C) to obtain 
a perturbation series for Z. Going over to the cumulants, 
we find the free energy per unit area 

/ = /m 2 + ^ J d 2 x(V int {(p{x))) m 2 c - 

/ d 2 xid 2 a; 2 (Mnt(^(xi))M n t(v3(x) 2 )) m 2 iC + . . . , 

(2.15) 
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where the subscript c indicate s the cum ulants . Inserting 
the expansion (^J) and using ( 2.10 ) and ( 2.12 ) , the series 
can be written as 



/ = m 
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(2.16) 



where the coefficients a n are dimensionless real numbers, 
starting with ao = 1/8. The higher expansion coefficients 
a n are combinations of integrals over the connected cor- 
relation functions: 



"2 = I rf2a; (</( x ))r 
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d 2 x 1 d 2 x 2 (<p 4 (xi)tp 4 {x 2 )) m 2 yC , 
m 12 /" 

e4£6 lT / d2xid2x2 (^ 6 ( X l) ^ 4 ( X 2))m 2 ,c 



(2.17) 



(2.18) 



(2.19) 



£4^4 J d 2 x 1 d 2 x 2 d 2 x 3 {ip A (x 1 )ip 4 (x 2 )(p 4 (x 3 )) m 2 c , 



To find the free energy ( 2.16 ) between walls, we must 
go to the limit m 2 — > 0. Following HQ, we substitute 
m 2 by the variational parameter M 2 , which is introduced 
via the trivial identity m 2 = y/ M 4 — <?r with r = (M 4 — 
m )/g, and expand this in powers of g up to the order 
g N . In the limit m 2 — ► 0, this expansion reads 

2/s ! r2\ 11^-2 1 r 1 7,2 2 1 7,3 3 

m [M ) = M r q ? o ^ q — ... . 

v y 2M 2S 8M« y 16M loy 

(2.20) 



Inserting this into (2.16), reexpanding in powers of g, 
and truncating after the iVth term, we arrive at the free 
energy per unit area 



f N (M 2 , d) = M 2 a b + Y, a n g n M 2 ^b n (2.21) 

n=l 

with 



N-n 



bn = £ (-1/ 



k ((l-n)/2 



k=0 



(2.22) 



being the binomial expansion of (1-1) (1 -")/ 2 truncated 
after the (N - n)th term Q. The optimization of Q2.21| ) 
is done as usual by determining the minimum of 
/n(M 2 , d) with respect to the variational parameter M 2 , 
i.e. by the condition 



df N (M 2 ,d) , 



9M 2 



= 0, 



(2.23) 



whose solution gives the op timal value Mjy(d). Resubsti 



tuting this result into Eq. ( [2.21 ) produces the optimized 
free energy /jv(<i) = fN(M^(d),d), which only depends 
on the distance as /iv(d) = ^a^jd 2 . Its derivative with 
respect to d yields the desired pressure law with the ,/Vth 
order approximation for the constant oln'- 



Pn = a N 



(2.24) 



We must now c alculate the cumulants occuring in the 
expansion (2.16). 



III. EVALUATION OF THE FLUCTUATION 
PRESSURE UP TO FOUR-LOOP ORDER 



The correlation functions appearing in (2.17)-( p.l9| ) 
are conveniently represented by Feynman graphs. Green 
functions are pictured as solid lines and local interactions 
as dots, whose coordinates are integrated over: 



xi 



x 2 = G m 2(xi,x 2 ), 



d 2 . 



(3.1) 
(3.2) 



These rules can be taken over to momentum space in the 
usual way. One easily verifies that t he in te grals over the 
connected correlation functions in ( 2.17 )— ( P ■ 1 9| ) have a 
dimension A/m 2 ^ n+l/_1 ^ , where V is the number of the 
vertices of the associated Feynman diagrams. Thus we 
parametrize each Feynman diagram by vA/m 2 ( n+v ~^ , 
with a dimensionless number v, which includes the mul- 
tiplicity. In Tab. |[ we have listed the values v for all 
diagrams up to four loops. No divergences are encoun- 
tered. Exact results are stated as fractional numbers. 
The other numbers are obtained by numerical integra- 
tion, which are reliable up to the last written digit. The 
right-hand column shows numbers vk obtained by the 
earlier approximation where all the Feynman dia- 
grams were estimated by an analogy to the the problem 
of a particle in a box. In Ref. B , it was shown that the 
value v of a large class of diagrams of the membrane prob- 
lem can be obtained by simply dividing the value of the 
corresponding particle-in-a-box-diagram vpb by a factor 
1/4 L , where L is the number of loops in the diagrams. 

Inserting the numbers in Tab. [j| into (2.17)— (2.19), we 
obtain the coefficients 01,02,(13 of the free energy per 
area (2.21), which is then extremized in M 2 . To see how 



the results evolve from order to order, we start with the 
first order 



h(M 2 ,d) = ^a Q M 2 



Oi 



d 2 



(3.3) 



with ao = 1/8 and a± = 1/64. Here, an optimal value of 
M 2 does not exist. Thus we simply use the perturbative 
result for m — which is equal to (3.3) for M = 0. Dif- 
ferentiating /i(0, d) with respect to d yields the pressure 
constant in (2.24): 
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0.038553. 



(3.4) 



This value is about half as big as the Monte-Carlo esti- 
mates ( |1.3| ) and agrees with the value found in M. To 
second order, the reexpansion (2.21) reads 



/ 2 (M 2 ,d) = ^a M 2 +oi J +a 2 ^^ (3.5) 

with a 2 « 1.0882 • 10~ 3 from Tab. | Minimizing this 
energy in M 2 yields an optimal value 



M 2 (d) 



8a 2 7T 2 
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0.152362 



d 2 ' 



(3.6) 



and 



Inserting ao = 1/8 and oi, a 2 from Tab. [|, we obtain 

a 2 w 0.073797, (3.8) 



thus improving drastically the first-order estimate ( |3.4| ) 
This value is by a factor 1.026 larger than that obtained 
in the approximation of Ref. || . 

Continuing this proceeding to third order, we must 



for the anharmonic oscillator fjl. Other systems treated 
by variational perturbation theory show a similar behav- 
ior 0. Assuming that an exponential convergence ex- 
ists also here, we may extrapolate the sequence of values 
ai,a 2 ,Q!3 calculated above to infinite order. It is useful 
to extend this sequence by one more value at the l owe r 
end, ao = 0, which follows from the one-loop energy (2.9) 
at rn 2 = 0. This sequence is now extrapolated towards 
a hypothetical exact value a ex by parametrizing the ap- 
proach as 



ajM = exp (— T) — £N e 



(4.1) 



The parameters 77, f, e, and the unknown value of a ex 
are determined from the four values ao, • • • , as, with the 
result 

r) = 2.529298, £ = 0.660946, e = 1.976207, (4.2) 



(3.7) and the extrapolated value for the exact constant: 



etc 



0.0797149. 



(4.3) 



This is now in perfect agreement with the Monte-Carlo 
values (|Q|). 

The approach is graphically shown in Fig. ^ where the 
optimized values ao, . . . , a 3 all lie on a straight line (solid 
line). For comparison, we have also extrapolated the first 
four values a^ , . . . , af< in the approach of Ref. B yielding 
a value a CX K ~ 0.0759786, which is 4.9% smaller than 



minimize 



(4.3) 



/ 3 (M 2 , d) = Aa M 2 + a, 1L + 3 fl2 ZL_L + a ^_L 

(3.9) 

with a 3 w 2.7631 • 10~ 5 . The optimal value of Af 2 is 



Mf(d) 



32 a 2 
5 a 



0.219608 
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(3.10) 



Inserted into (3.9), we find the four-loop approximation 
for the proportionality constant a: 
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a ex - an 



a 3 w 0.079472. 



(3.11) 



This result is i n very good agreement of the Monte-Carlo 
results in (1.3). It differs from the approximate value of 
the method presented in Ref. Qj by a factor 1.047. 

An even better result will now be obtained by extrap- 
olating the sequence a%, a 2 , a 3 to infinite order. 



FIG. 1. Difference between the extrapolated pressure 
constant cv ox and the optimized N-th order value qjv ob- 
tained from variational perturbation theory for the method 
presented in this paper (solid line) and the first four values 
of the approximation scheme introduced in Ref. [^] (dashed 
line). Dots represent the values to order N in these approxi- 
mations. 



IV. EXTRAPOLATION TOWARDS THE EXACT 
CONSTANT 

Variational perturbation theory exhibits typically an 
exponentially fast convergence. This was exactly proven 



V. SUMMARY 



We have calculated the universal constant a occur- 



ing in the pressure law (1.1) of a membrane fluctuat- 
ing between two walls. This has been done by replacing 



4 



the walls by a smooth potential with a parameter m 2 . 
This potential approaches the wall potential in the limit 
m 2 — > 0. The anharmonic part of the smooth poten- 
tial was treated perturbatively. The limit m 2 —> cor- 
responds to a strong-coupling limit of the power series, 
and was calculated by variational perturbation theory. 
Extrapolating the lowest four approximations to infin- 
ity yileds a pressutre constant a, which is in very good 
agreement with Monte-Carlo values. 
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TABLE I. Feynman diagrams with loops L, multiplicities 
s, and their dimensionless values v. The last column shows 
the values i>k = i>pb/4 l used in Ref. M. 
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